Landau-Zener-Stiickelberg interference in a multi-anticrossing system 
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We propose a universal analytical method to study the dynamics of a multi-anticrossing system subject to 
driving by one single large-amplitude triangle pulse, within its time scales smaller than the dephasing time. Our 
approach can explain the main features of the Landau-Zener-Stiickelberg interference patterns recently observed 
in a tripartite system [Nature Communications 1:51 (2010)]. In particular, we focus on the effects of the size of 
anticrossings on interference and compare the calculated interference patterns with numerical simulations. In 
addition, Fourier transform of the patterns can extract information on the energy level spectrum. 

PACS numbers: 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

Landau-Zener (LZ) transitions^ at the anticrossings play 
a fundamental role in coherent quantum control which is key 
to the realization of quantum computation 4 . To date, great ef- 
forts have been devoted to investigate the coherent quantum 
dynamics of the states at energy-level crossings. In a strongly 
harmonic-driven two-level system (TLS), repeated LZ transi- 
tions give rise to Stiickelberg or Ramsey-type oscillations, in 
analogy to Mach-Zehnder (MZ) interferometer. The MZ-type 
interferometry has been observed in a driven superconducting 
flux qubit^, a Cooper-Pair box-, and a quantum dot system 7 . 
These patterns have been theoretically reconstructed from dif- 
ferent perspectives^—. 

In addition to Landau-Zener-Stiickelberg (LZS) interfer- 
ence of one-anticrossing type, several experiments have also 
reported LZS interference in the multi-anticrossing level 
structure. In the presence of much stronger harmonic exci- 
tation, the qubit state can be driven through more of the con- 
stituent energy-level anticrossings, and the resulting LZS in- 
terference reveals complex checkerboardlike patterns^. This 
has been explained well in theor y 13 ' 14 . Most recently, a new 
method of coherent manipulation of quantum states in a tri- 
partite quantum system formed by a superconducting qubit 
coupled to two TLSs is reported^. The manipulation relies 
on the LZS interference produced by transitions at the two an- 
ticrossings, and has potential application in precise control of 
quantum states in the tripartite system. Nevertheless, a univer- 
sal model to explain the observation has not been proposed. 
We shall show below that it can be understood rather easily 
using our approach. 

In this work, we start with a strongly driven one- 
anticrossing system where we consider LZ transition as a 
gate operation. We then analyze the dynamics of multi- 
anticrossing system. Under the strong triangle pulse driv- 
ing, occupation probability of the system at the initial state 
exhibits, as a function of the driving amplitude and the 
pulse width, diverse LZS interference patterns. Our approach 
presents a unified analytical treatment of these diverse pat- 
terns. In a specific case of two-anticrossing system, we focus 
on how the interference patterns are influenced by the sweep 



rate and the coupling strength, and thereby elucidate the un- 
derlying physics of the various patterns. Converting the pat- 
terns to the phase domain, Fourier transform of the resulting 
population oscillation reveals Fourier components of the com- 
pound pattern, which is in agreement with our analysis. In all 
the systems under study, the influence of relaxation and de- 
phasing is neglected to obtain a clear physics picture of the 
underlying quantum physics. 

This paper is organized as follows. In Sec. II we intro- 
duce the analytical expression to describe the dynamics for 
the multi-anticrossing system. In Sec. Ill we apply the general 
result to the N = 2 case, and focus on the effect of the size 
of anticrossings on the interference patterns. The calculated 
patterns of four representative combinations of two anticross- 
ings are in agreement with the numerical results. A discussion 
of the formation of the interesting dark state-one special case 
in the two-anticrossing system, is also included. In Sec. IV, 
Fourier transform of these patterns are presented, exhibiting 
an explicitly ordered structure of one dimensional arcs that 
offers energy spectrum information. Finally, Sec.V contains a 
summary of this work. 



II. ANALYTICAL MODEL 

First, we consider one-anticrossing system in which two 
energy levels of a quantum TLS "cross" each other as some 
external parameter is varied. At the energy-level crossing, hy- 
bridization of the two states results in an anticrossing due to 
coupling of the states, as shown in FigJTjA). The Hamiltonian 
of TLS is, 



Htls 



e (t) A 
A ei (f) 



(1) 



where erj(i) and ei(i) are the energy levels of two diabatic 
states, and A is their coupling strength. The transition be- 
tween energy levels at the anticrossing is what we call LZ 
transition. 

Following Damski and Zurek's adiabatic-impulse approxi- 
mation model 3 -^, we can obtain a convenient description of 
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FIG. 1: (Color on line) Schematic diagram of (A) one anticrossing 
and (B) LZS interference in a phase qubit coupled to N TLSs. (A) 
LZ transition at the anticrossing splits the initial state into a super- 
position of two states in a like manner of a beam splitter. (B) Start- 
ing at the dot maker, the qubit state is swept by a triangle pulse. A 
succession of LZ transitions at the chain of anticrossings result in a 
superposition state. The state accumulates phases (shaded region), 
interferes at the return LZ transitions successively, and returns to the 
initial state with a finite probability. 



the system's dynamics. It is provided that the system evolves 
adiabatically everywhere except at the points of minimum en- 
ergy splitting where a sudden mixing in the population of the 
two energy levels occurs. This non-adiabatic transition at the 
anticrossing can be described by a unitary transformation^ . 



cos(0/2) exp(-i<ps) isin(0/2) 

isin(0/2) cos(0/2) exp(i</> s ) 



(2) 



where sin 2 (0/2) = Plz, where Plz is the Landau-Zener 
transition probability at the anticrossing. If the anticrossing is 
swept from the infinity on one side to the infinity on the other, 
Plz has the asymptotic form: 



A 2 

Plz = cxp(-2tt— ), 
nv 



(3) 



in which v = d{e\ — eo)/dtis the variation rate of the energy 
separation between the two diabatic levels, and 2A is the size 
of the anticrossing. In addition, phase jump 4>s = (j>s ~ 7r /2, 
related to the general Stokes phenomenon^. <f>s is the so-called 
Stokes phase that takes the form: 



where 8 = A 2 /his is called adiabatic parameter and T is the 
gamma function. In the adiabatic limit <fis — > 0, and in the 
sudden limit <f>s — >• 7r/4. 

From the perspective of optics, the avoided level crossing, 
when driven through, can be viewed as a beam splitter, be- 
cause LZ transition taking place at the anticrossing splits an 
input state in a superposition of two states, just analogous to 
an optical beam splitter which splits the incident light in two. 
In this sense, we define reflection coefficient \r\ 2 = 1 — Plz 
and transmission coefficient \t\ 2 = Plz- 

Then, we take into account the multi-anticrossing system 
which can be realized in a superconducting phase qubit with 
many TLSs inside its Josephson Junction. A phase qubiti^r— 
consists of a single current-biased Josephson junction. When 
biased close to the critical current Iq, the qubit can be treated 
as a tunable artificial atom with discrete energy levels that ex- 
ist in a potential energy landscape determined by the circuit 
design parameters and bias. In the qubit-TLS coupled system, 
TLSH~— is formed in the disordered barrier material, where 
some atoms can occupy two positions, corresponding to two 
quantum states^. When the energy separation of qubit states 
equals that of TLS, resonant tunneling between the states 
opens an avoided level crossing in the energy spectrum of 
qubit, and forms a multi-anticrossing chain (see FigQlB)). Its 
Hamiltonian takes the form, 
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(5) 



The time dependent e(t) is energy spacing between the ground 
state and the excited state of qubit and i = 1, • ■ • , N is en- 
ergy spacing of the i th TLS. It is assumed that qubit is initially 
prepared in its excited state, and quantum state transitions are 
then driven using a triangle pulse with amplitude V and pe- 
riod T. This is a double-passage process during which anti- 
crossing regions are passed twice. The first excursion through 
the anticrossings coherently divides the signal into N output 
paths, where the dynamical phase is accumulated during the 
adiabatic parts of the evolution. Then the second excursion re- 
combines the separated signals via LZ transition, and results 
in the LZS interference patterns. 

The algebra describing this one full driving cycle of triangle 
pulse is as follows. As the anticrossings are traversed from left 
to right in turn, according to Eq.© the transition amplitude to 
each output path is: 

' cos 0i to path 1 

sin 0i cos 02 to path 2 

sin 0i sin 2 cos 03 to path 3 



-|-5(ln5- l) + axgT(l-iS), 



(4) 



(nit 



(ntr 



sin 9k ) cos 8i to path i 



(6) 



j. ri, — _l 

nk=N . 
fc=l Sm 



sin 0fc ) cos 0at to path N 

to path N + 1 



3 




FIG. 2: (Color online) Different combinations of two beam-splitters. (A) A a < As- Mi acts as a beam splitter, while M2 reflects the state 
with Plz approaching zero. Thus the interfering pattern is mainly dependent on the phase accumulated in region I. (B)We still set A a < As, 
but both of them are much smaller so that Mi transmits the state with Plz approaching unity, while M2 plays the role of a beam splitter. 
The phase accumulated in region II plays a major role in forming the interference fringes. (C)Aa > As- Mi acts as a beam splitter, while 
M2 transmits the state. The phase accumulated in region I+II plays the leading part. (D)Aa — As, Both Mi and M2 act as beam splitters. 
"3?/, ($/ + altogether contribute to the interference pattern. 



The adiabatic evolution is in the form of exp(i<pi), where </>j 
denotes the phase accumulated on path % relative to the ground 
state |0) between the two successive crossings. When the anti- 
crossings are traversed from right to left, the transition ampli- 
tude back to the initial state from each path Ai n equals A out . 
Then after one round of LZS interference, the amplitude for 
the coupling system remaining in the initial state is: 

pathN+l 

Atotai = 2J A out (i)Ai n (i) exp(i(/>i) 

i—pathl 
pathN+l 

= A LWexp(#i). 

i—pathl 

Therefore, the probability for the upper level after one period 
of triangle pulse is given by: 

Pi = \A tota i\ 2 

N+l N+l 

= A lut{i)A 2 ut{3)eMi4>ii)- 
i=i j=i 

It can be further simplified to the form, 

i=N+l 

Pi= E A ™&) 

i=l 

7V+1 j<i 

»=i j=i 

where 

n—i — 1 

faj = 4>i - h = E ®n-(j<i) (7) 

n=j—l 

Here $„ = J^[E pathn (t) - E pathn+1 (t)]dt (at t = t n (t' n ), 
the 77th anticrossing is traversed from left (right) to right (left)) 



is the interference phase accumulated in the area of region n 
as shown in Fig[TjB). 

Eq.((7]i indicates that the resultant interference pattern is 
subject to the LZ transition amplitude and the interference 
phase. Every two paths in the energy diagram accumulate 
one interference phase and give rise to one type of interfer- 
ence fringes (see Fig[TjB)). The weightings of the C^r +1 in- 
terference patterns governed by a single phase difference on 
the resultant interference depend on LZ transition amplitude. 
The two factors to determine the LZ transition probability 
are the size of the anticrossing and the velocity with which 
it is traversed. Therefore by varying sweep rate and coupling 
strength between qubit and TLSs, we can generate a variety of 
interference patterns with promising use in quantum control. 

For example, (1) sin#fc ~ 0, which means the size of the 
kth anticrossing is so large that no ingredient of the wave- 
function can transmit through it. Therefore, for i > k, 
A ou t{i) = 0, and the number of interferences reduces from 
Cjj +1 to Cf . Especially, if k = 2, there is totally Cf = 1 
interference, i.e., the interference between path 1 and path 2. 
(2) cos 8k — 0, which means the size of the kth anticrossing 
is so small that the wavefunction cannot feel its existence. In 
this case, A out (k) = 0, and the kth path will not participate 
in the interference. For more interesting concrete examples, 
let's turn to Sec. III. 



III. APPLICATION TO TWO ANTICROSSINGS 

LZS interference in a phase qubit coupled to two TLSs has 
recently been observed^. 3D view of the probability for the 
initial state of qubit under one round of strong triangle pulse 
drive explicitly characterizes the sweep-rate-dependency of 
LZS interference. This feature can be explained using the 
above result applied to N=2 case. In this case, the coupled 
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FIG. 3: (Color on line) (A)-(D) are numerical simulations of qubit population at the initial state|l), plotted as a function of the amplitude of the 
driving triangle pulse A pu i se and its time width T pu ise- (E)-(H) are analytical results using Eq.([9](. The parameters used are extracted from our 
pertinent experiment. (A & E)A A = WMHz,A B = 60MHz. (B & F)A A = 1MHz, A B = WMHz. (C & G)A A = WMHz, A B = 
1MHz. (D & H)Aa = As = 17MHz. 




FIG. 4: (Color on line) Discrete Fourier Transform of the LZS patterns in Fig[3] In FT, one curve corresponds to one interference pattern. The 
number of the curves indicate the number of anticrossings that act as beam splitters under the pulse driving. The curves in (A) (B) (C) indicate 
<&n, &i + <&n pattern respectively. (D) shows all the three curves as in this case all the patterns make comparable contributions to the 
interference. 



Hamiltonian is, 

Hqubit-2TLSs = 



e(t) Ax A 2 
Ai ei 
A 2 e 2 



(8) 



According to the general formula in Eq.©, the occupation 
probability at the initial state Il.91.g2) after one pulse driving 
takes the form: 



sin 4 9\ cos 4 62 



sin 4 61 sin 4 9 2 



Pi = cos 4 01 

+ 2 sin 2 61 cos 2 6 2 cos 2 9 l cos <£/ 

+ 2 sin 4 6»i sin 2 6» 2 cos 2 6 2 cos $7/ 

+ 2 sin 2 6»i sin 2 2 cos 2 B\ cos($ 7 + $//), 



(9) 



where <&;(i=I,II) is the total phase accumulated in region i as 
shown in Fig. 2. It clearly exhibits that the interference fringes 
comprise of 3 (Cf) patterns, governed by phase accumulation 
in area I, II, and I+II, respectively. With the size of the two 
anticrossings fixed, as the sweep rate is varied, the interfer- 
ence pattern can be divided into three regions displaying the 
features of the three phase patterns respectively. In the slow 
limit, TLSi acts as a beam splitter and TLS 2 acts as a total re- 
flection mirror; $/ pattern dominates. In the fast limit, TLSi 
acts as a total transmission mirror and TLS 2 acts as a beam 
splitter; pattern makes the main contribution. In the in- 
termediate region, both TLSi and TLS 2 act as beam splitters; 
<f>j + <J>/j pattern is the principal feature. 

On the other hand, the size of anticrossing also manipulates 
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weights of the three patterns in the resulting compound inter- 
ference fringes. Though experimental realization of coupling 
strength control has not been achieved, investigation on the 
effect of the size of anticrossings on the interference patterns 
could make predictions of various interference fringes with 
potential use in future coherent control of hybrid qubit sys- 
tem. Here we choose four representative combinations of two 
TLSs, and the parameters used are extracted from G.Z.Sun's 
experiment 1 -. Based on this we discuss the effect of A: 

(i) A A = 10MHz, A B = 100MHz. In this case, P LZ2 ~ 
0, sin #2 ~ 0, cos 02 ~ 1. The occupation probability approx- 
imately equals to: 




Pi = sin 4 6»i 



cos 



6»i +2 sin 2 6»i cos 2 6»i cos $/. (10) 



FIG. 5: Dark state in a hybrid qubit. (A) The two TLSs share the 
same location in the energy diagram. (B) Calculated energy spec- 
trum from Hamiltonian d 13t.The red branch keeps staying at the state 
\TLS1). 



This combination reveals the main feature of $ j pattern. 

(ii) A A = lMHz,A B = 10MHz. In this case, P LZ1 ~ 
1, sin 6*i ~ 1, cos6>i ~ 0. The occupation probability approx- 
imately equals to: 



Pi 



sin 4 2 + cos 4 02 + 2 sin 2 2 cos 2 2 cos $// . (11) 



This combination reveals the main feature of $// pattern. 

(iii) A A = 10MHz, A B = 1MHz. In this case, P LZ2 ~ 
1, sin 2 ~ 1, cos 02 ~ 0. The occupation probability approx- 
imately equals to: 



Pi = sin 4 0i + cos 4 0i + 

+ 2 sin 2 0i cos 2 0i cos($/ + $//) 



(12) 



This combination reveals the main feature of pattern. 

(iv) Aa = A B — 10MHz. The occupation probability is 
in the form of Eq.©. 

Provided that LZ transition probability takes the asymptotic 
form in Eq.([3j, interference patterns in the 4 cases above cal- 
culated based on Eq.® (see Figj3|A)-(D)), agree well with 
the numerical simulations (see Figf3tE)-(H)), except the slight 
modulation in one interference fringe involved in the numeri- 
cal results. The modulation is caused by the fluctuation of ac- 
tual LZ transition probabilities around its asymptotic form 1 ^. 

In addition, it is noteworthy that in a qubit-two-TLSs hy- 
brid system, it is completely possible there exists a smaller 
anticrossing screened by a larger one. Although the LZS in- 
terference is just the same as that of a single anticrossing, one 
branch in the spectrum is always at the excited state of TLS, 
\TLS1) (see red line in Fig|5lB)). The system's simplified 
Hamiltonian in a basis formed by \lgig 2 ), [Oei^), |0<?xe2) 
reads, 



Ht 



w Oi/2 2 /2 
fii/2 

n 2 /2 



(13) 



where u> is the detuning between qubit and TLS. The eigen- 
state corresponding to the particular branch is: 



\*d) 



|0eiff 2 ) 



|0 5l e 2 ). (14) 



With no ingredient of 1 19132) in it, the state \$>d) is analogous 
to the Dark State in quantum optics, which has no ingredient 



of excited state. In this sense, we can call this "Dark state of a 
hybrid qubit". Vacancy of |l<7i<?2) state could largely reduce 
the influence of environment on the dark state, and thereby 
has potential use in information storage. The specific method 
to realize this is in need of further investigation. 



IV. FOURIER TRANSFORM 

Furthermore, Fourier transform is a helpful tool in identi- 
fying the individual phase components of the compound LZS 
interference. Supposing the two-anticrossing system follows 
a linear ramp traversing the anticrossing regions, FT of the 
occupation probability in Eq.© is found to be, 

P FT (kT,A pu [ se ) = B + B\5(k T - k T i{A pulse )) 
+ B 2 S(k T - k T2 (A pu i se )) 
+ B a S(k T - k T i(A pu [ se ) - k T2 {A pulse )), 

(15) 



where 



Bq = cos 4 0i + sin 4 0i cos 4 02 + sin 4 0i sin 4 02 
Pi = 2 sin 2 0i cos 2 2 cos 2 0i, 
P 2 = 2 sin 4 0i sin 2 2 cos 2 2 , 
P 3 = 2 sin 2 0i sin 2 2 cos 2 0i, 



kri = £12 
k T2 



(ex + e 2 )ex 



*2sA pu i se 
{sA pu i se — ei 2 ) 2 
2sA pu i se 



Here s is the diabatic energy-level slope, and £12 is the detun- 
ing between the location of anticrossing I and II. 

Therefore, it is expected that FT reveals a highly ordered 
structure of one-dimensional arcs in Fourier space. (A) 
pattern dominates, i.e., sin 02 ~ 0. In this case, P2, P3 ~ 0, 
and Pft — 2 sin 2 81 cos 2 8iS(kT — kTi(A pu i se )). Appar- 
ently, there is only one curve corresponding to kp = fcrx 
after FT. (B) $jj pattern dominates, i.e., cos0i ~ 0. In this 
case, Px, P3 — 0, and Pft — 2 sin 4 0x sin 2 02 cos 2 9 2 8(kT — 
kT2(A pu ise))- Only kr = &T2 shows up after FT. (C) 
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$j + <5ji pattern dominates, i.e., cos 62 — 0. In this 
case, B 1 ,B 2 ~ 0, and P FT ~ 2 sin 2 ^ cos 2 9 1 S(k T — 
kTi{A pu i se ) - k T2 {A pu i se )). Only k T = k T3 shows up after 
FT. (D) All three patterns including $/, $// and $/ + $// 
dominate, i.e., sin 2 #1 ~ sin 2 62 — 1/2. In this case, 
B\ ~ i? 3 ~ 1/4 and B 2 1/8. All three curves corre- 
sponding to kx = hri, kr = &T2 and kx = &T3 can be 
observed after FT. 

This is explicitly demonstrated in Figj4] which is based 
upon the Discrete Fourier Transform method: 

N 

P D F T (k T ,A pulse ) ^Y.P^^p^)^ 1 ^ 1 ^ < 16 > 
i=i 

where cjjv = exp(—2iri/N),N = 1,2,.... As labeled in the 
figure, different curves correspond to different phase patterns, 
which is in good agreement with our analysis above. In the 
multi-anticrossing system, if N anticrossings take part in the 
LZS interference, C^ +1 phase components can be observed 
in its FT and vice versa. Therefore, FT provides a means to 
ascertain how many TLSs are effectively coupled to qubit. 

V. SUMMARY 

In summary, we present a simple form of analytic expres- 
sion for describing the coherent dynamics of a driven multi- 
anticrossing chain. The oscillatory population of the hybrid 
system remained at the initial state exhibits a rich pattern of 
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